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ABSTRACT 

We analyze a sample of 82850 stars from the RAVE survey, with well-determined 
velocities and stellar parameters, to isolate a sample of 18026 high-probability thin-disc 
dwarfs within 600 pc of the Sun. We derive space motions for these stars, and deduce 
the solar space velocity with respect to the Local Standard of Rest. The peculiar solar 
motion we derive is in excellent agreement in radial Uq and vertical Wq peculiar 
motions with other recent determinations. Our derived tangential peculiar velocity, 
Vq agrees with very recent determinations, which favour values near 13 km s -1 , in 
v> , disagreement with earlier studies. The derived values are not significantly dependent 

j_j ■ on the comparison sample chosen, or on the method of analysis. The local galaxy 

seems very well dynamically relaxed, in a near symmetric potential. 

Key words: Galaxy: kinematics and dynamics, solar neighbourhood, solar motion, 
local standard of rest, stellar content 



1 INTRODUCTION 

The kinematics of stars near the Sun provide vital informa- 
tion regarding the structure and the evolution of the Galaxy. 
Since the observations are relative to the Sun, which is not 
stationary, the kinematical parameters we obtain need to be 
corrected accordingly to the Local Standard of Rest (LSR), 

* E-mail: basarc@istanbul.edu.tr 



which is defined as the rest frame of a star at the location of 
the Sun that would be on a circular orbit in the gravitational 
potential of the Galaxy. This implies the determination of 
the solar peculiar velocity components relative to that stan- 
dard, Uq, Vq and Wq. These values are also of considerable 
intrinsic interest, as their determination as a function of the 
spatial distribution and kinematic properties of the com- 
parison stellar sample is a test of the symmetry of the local 
Galactic potential as a function of, for example, velocity dis- 
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persion. One can test if thin disk slowly moving stars are, or 
are not, affected by a Galactic bar, among many other con- 
siderations. Uncertainty in the solar kinematic parameters 
can have pr ofound implications for Gal actic structural anal- 
yses (cf. e.g. lMcMillan fc Binnev"ll2010l ) and may be affected 
systematically in turn by presumpt ions about larger scal e 
Galactic stellar populations (cf. e.g. ISchonrich et akll201Cj ). 
There have been several very recent determinations of the 
solar peculiar velocity, including some using RAVE (RAdial 
Velocity Experiment, Stein metz et al., 2006, see also Section 
2) data as we use here (e.g. IPasetto et alj|201(ih . A summa- 
rized list of recent determinations of the space velocity com- 
ponents of the Sun appearing in the lite rature is given in 
Table 1, while iFrancis fc Anderson] (|2009l ) provide an even 
longer list of references. Given the astrophysical significance 
of systematic sample- and method-dependent effects in the 
determination, we consider here a further determination of 
the solar space velocity, using a complementary technique. 

The determination of the solar space motion relative 
to a Local Standard of Rest (LSR) is a long-standing chal- 
lenge. The robust determination of the systematic motion of 
some stellar population relative to some agreed LSR would 
itself be of considerable intere st. Hence many studies have 
been carried out since that of iHomannl l| 18861) to estimate 
these parameters. However, there is still uncertainty about 
the numerical values of the solar velocity components, rela- 
tive to local thin disc stars, especially for the Vq component, 
due to the complexity involved in compensating for the ve- 
locity "lag" of sample stars. The difference between the re- 
sults of different researchers originates from different data as 
well as diffe rent procedures us e d. Th e very different results 
derive d by iDehnen fc Binnevl (|l998l ) and ISchonrich et al.l 
(|2010h . both of whom who used the Hipparcos data, are 
an illustration of the continuing model-dependence of the 
result. 

The radial (Uq) and vertical (Wq) components of the 
solar motion, as a function of the relevant stellar popula- 
tion, can be obtained in principle by direct analysis of the 
mean heliocentric velocities of solar neighbourhood stars, 
subject to appropriate sample selection. As one can see in 
Table 1, the various determinations are in tolerable agree- 
ment. However, the determination of the component in the 
direction of Galactic rotation (Vq) is complicated by sys- 
tematics: the mean lag, i.e. asymmetric drift V a , with re- 
spect to the LSR which depends on the velocity disper- 
sion (a) of the stellar sample in question. An extensive 
discussion of the radial Jeans' equation on which analysis 
of the asymmetric drift is based, including its derivation, 
relevant approximations, and app lications, can be found in 
iGilmore. Wvse fc Kuiikenl (|l98Sl S- At its simplest, a single 
population asymmetric drift/ Jeans' analysis assumes a lin- 
ear relation between (suitable values of) the asymmetric 
drift V a of any stel lar sample and it s squared radial veloc- 
ity dispersion o\ (|Stromberg| 1 19461 ). Formal least-squares 
straight line fits intercept the V a axis at some negative 
value, typically in the range -20 to -5 km s -1 . This V a in- 
tercept determines the solar velocity Vq using the equation 
Vs=V a +VQ. Vs is the negative mean heliocentric azimuthal 
veloci ty of any stellar sample. T he value Vq=5.25±0.62 km 
s _1 of lPehnen fc Binnevl (|l998l ) (hereafter DB98) is typical 
of the smaller values deduced recently from this procedure. 

However, four other recent studies argue for larger val- 



ues of Vq than that of DB98. IPiskunov et""ai1 (|2006l ) de- 
termined a value different from that of DB98 by about 7 
km s _1 in the V component of the solar motion, when de- 
termined with respect t o open clusters in the solar neigh- 
bourhood. iBinnevI (|2010| ) fitted distribution function models 
to two sets of velocity distributions and obtained a differ- 
ence of about 6 km s _1 compared to the valu e of DB98. 
This higher value of Vq has been confirmed by Rcid et al.l 
(|2009h and iRygl et all (|2010h in their works related to ra- 
dio frequency astrometry of masers in regions o f mas sive 
star formation (thou gh see McMillan fc Binnevl (|20 10T ) for 
a different analysis). ISchonrich et al.l (|2010l ) applied a par- 
ticular chemodynamical model of the Galaxy to DB98's data 
and determined a value for the Vq solar velocity component 
consistent with recent high values. The key factor in the 
ISchonrich et~ail (|2010l ) analysis, and of others, is the rele- 
vant range of radial v elocity dispersions cr 2 ? which is used 
in the Jeans' analysis. ISchonrich et al.l (|2010l ) showed that 
for crf(^600 (km s -1 ) 2 the Hipparcos data define a straight 
line, whereas for erf^400 (km s -1 ) 2 they deviate from such 
a simple fit. If one omits stars with low dispersions, and ex- 
trapolates the linear fit for stars with high dispersion to zero, 
the linear fit for stars with high dispersions, one confirms the 
low value of DB98. Actually, this is the procedure used by 
DB98, who also ignored stars with low velocity dispersion 
due to their probable lack of dynamical equilibrium. While 
it is true that low velocity dispersion stars are likely to be 
most affected by the effects of dissolving star clusters and 
the non-axisymmetric gravitational potential of spiral arms, 
one does need to consider carefully the offset character of 
the relation between V a and o\ at low velocity dispersions 
to avoid bias. 

Here we use a new sample of stars and a different 
methodology to re-estimate the Uq, Vq and Wq solar veloc- 
ity components. 1) We use RAVE data extending to larger 
distances than the Hipparcos sample. 2) We applied the fol- 
lowing constraints to obtain a sample of main sequence stars: 
i) we selected stars with surface gravity 4<log<?<5, ii) we 
omitted stars with (J — H) <0.05 and (J — H) >0.55 to 
avoid the blue horizontal branch and possible red giant stars, 
iii) we excluded stars with space velocity errors larger than 
25 km s _1 , iv) we separated stars into populations (see Sec- 
tion 5) and we used the thin disc population which is not 
contaminated by thick disc/halo as a preferred sample in 
our work. The last constraint is especially important in es- 
timating Vq, because it excludes thick disc and halo stars 
which have relatively large asymmetric drift. Thus, it limits 
the range of V velocity component, which in turn minimizes 
possible astrophysical large-scale dynamical asymmetry ef- 
fects. 

The RAVE survey is described in Section 2, the data 
are presented in Section 3, Sections 4 and 5 are devoted to 
kinematics and population analysis, respectively. The results 
are given in Section 6 and finally a discussion is presented 
in Section 7. 



2 RAVE 



The RAdial Velocity Experiment (RAVE, ISteinmetz et al.l 
2006) is a spectroscopic survey aiming to measure radial 
velocities and stellar atmospheric parameters, temperature, 



© RAS, MNRAS 000, [[]-?? 



LSR from RAVE DR3 3 



Table 1. Space velocity components of the Sun with respect to the LSR as given in the literature. 



Reference 


Source 


Uq 




Wq 






(km s 


(km s *) 


(km s — ) 


This study (2010) 


RAVE DR3 


8.50±0.29 


13.38±0.43 


6.49±0.26 


Bobylev & Bajkova (2010) 


Masers 


5.5±2.2 


11. Oil. 7 


8.5±1.2 


"r>„ u_i_ j ,1 Com r\1 

Breddels et al. 


HAVE DK2 


12.0±0.o 


20.4±0.5 


7.8±0.3 


bchonrich et al. (2010) 


Hipparcos 


11.10±0.72 


12.24±0.47 


7.25±0.36 


Francis k, Anderson (2009) 


Hipparcos 


7.5±1.0 


13.5±0.3 


6.8±0.1 


Veltz et al. (2008) 


RAVE DR1 


8.5±0.3 




ll.lil.O 


Bobylev & Bajkova (2007) 


F & G dwarfs 


8.7±0.5 


6.2±2.2 


7.2±0.8 


Piskunov et al. (2006) 


Open clusters 


9.44±1.14 


11.90±0.72 


7.20±0.42 


Mignard (2000) 


K0-K5 


9.88 


14.19 


7.76 


Dehnen & Binney (1998) 


Hipparcos, d ma x = W0 pc 


10.00±0.36 


5.25±0.62 


7.17±0.38 


Binney et al. (1997) 


Stars near South Celestial Pole 


11±0.6 


5.3±1.7 


7.0±0.6 


Mihalas k Binney (1981) 


Galactic Astronomy 2 nd Ed. 


9.2±0.3 


12.0 


6.9±0.2 


Homann (1886) 


Solar neighbourhood stars 


17.4±11.2 


16.9±10.9 


3.6±2.3 



metallicity, surface gravity, of up to one million stars using 
the 6 degree Field (6dF) multi-object spectrograph on the 
1.2 m UK Schmidt Telescope of the Australian Astronomi- 
cal Observatory. The RAVE programme started in 2003, ob- 
taining medium resolution spectra in the Ca-triplet region 
(8410-8795 A) for southern hemisphere stars in the mag- 
nitude range 9<Idenis<13. The scientific goals of RAVE 
include analyzing the chemical and dynamical evolution of 
the Galaxy, using both dwarfs and giants observed locally. 
The main sequence stars occupy a region extending to a dis- 
tance of a few hundred parsecs, whereas giants extend up 
to a few kpc. RAVE was designed to avoid high reddening 
in the Galactic plane and contamination from the bulge, i.e. 
|fo|>5° and £<315° constraints are applied when selecting 
programme stars. 

RAVE is a precursor of the spectroscopic part of the cor- 
nerstone mission Gaia of the European Space Agency. The 
wavelength range for RAVE spectra was cho sen to match 
that of the Gaia Radial Velocity Spectro meter (|Munaril2003l ; 
iKatz et al. I l2004l ; I Wilkinson et al.l |200i), i.e. around the Ca 
II IR triplet. This wavelength range also includes lines from 
elements such as Fe, Ca, Si, Mg and Ti which can be used 
to estimate [q/Fe] in addition to overall metallicity (see 
ISteinmetz et alj|2006l . for a more detailed description of the 
goals of RAVE). 



3 DATA 

The data used in this study are a working version of what 
will become RAVE's third data release (DR3; Siebert et al., 
in prep.). DR3 will consist of 82850 stars, each with equa- 
torial and Galactic coordinates, radial velocity, metallicity, 
surface temperature and surface gr avity. We also note th e 
two e xisting data re l eases , i.e. DR1 |Steinmetz et al.l [2006). 
DR2 ( Zw itter et al.l l200St ) . Proper motions were compiled 
from several catalogues: Tycho-2, Supercosmos Sky Survey, 
PPMXL and UCAC-2. Proper motion accuracy decreases in 
this order, therefore, if proper motions were available from 
all catalogues, Tycho-2's value was used. If Tycho-2 did not 
provide proper motions the values were taken from Super- 
cosmos Sky Survey, etc. Photometric data are based on op- 
tical and near infrared (NIR) systems. The magnitudes of 



stars were obtained by matching RAVE DR3 with Tycho-2, 
USNO-B, DENIS and 2MASS catalogues. The analysis here 
uses stars in RAVE DR3 with 2MASS catalogue photome- 
try. 



3.1 Main Sequence Stars and Distance 
Determination 

We applied two constraints to obtain a main sequence sam- 
ple: we selected stars with surface gravities 4<log<?<5 and 
we excluded stars with (J — i/)o<0.05 and (J — iif)o>0.55. 
The second constraint is especially effective in reducing the 
contamination due to blue horizontal branch and red gi- 
ant stars. Thus, the sample was reduced to 21310 stars. The 
(J — H)o-{H — K s )o two colour diagram of the parent sample 
(Fig. la) has a bimodal distribution, confirming the presence 
of K and M giants in addition to the F-M dwarfs, whereas 
the sample obtained after applying the two constraints cited, 
shown in Fig. lb, has only one mode indicating a pure main 
sequence sample. 

Contrary to the Hipparcos catalogue, parallaxes are not 
available for stars observed in the RAVE survey, but not ex- 
isting in the Hipparcos survey. Hence, the distances of the 
main sequence stars were calculated using another proce- 
dure. We applied t he ma in-sequence colour-luminosity rela- 
tion of iBilir et all {2008), which, as described in that refer- 
ence, is valid in the absolute magnitude range 0<Mj<6. The 
errors of the distances were estimated combining the inter- 
nal errors of the coefficients of Bilir et al.'s (2008) equation 
and the errors of the 2MASS colour indices. 

As most of the stars in the sample are at distances larger 
than 100 pc, their colours and magnitudes will be affected 
by interstellar reddening. Hence, distance determination is 
carried out simultaneously with de-reddening of the sample 
stars. As the first step in an iterative process, we assume 
the original (J — H) and (H — K s ) colour indices are de- 
reddened, and evaluate the Mj absolute magnitudes of the 
sample sta r s by means of the colour-luminosity relation of 
I Bilir et all (120081 ). Combination of the apparent and abso- 
lute magnitudes for th e J band gives the distance d of a star . 
We used the maps of ISchlegel, Finkbeiner fe Davis! (|l998l ) 
and evaluated the colour excess Eoo{B — V) for each sample 
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Figure 2. E(B — V) reddening map in the solar neighbourhood for our RAVE main sequence sample, (a) stars with \b\ ^ 30°, (b) stars 
with 30° < |6| < 60° and (c) stars with |6| > 60°. 
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Figure 1. Two colour diagrams for two samples of RAVE DR3 
stars, (a) for the original sample (main sequence stars and giants), 
(b) for main sequence stars only. 



star. The relation between the total and selective absorp- 
tions in the UBV system, i.e. 



A x (b) = 3AE X (B -V) 



(1) 



gives Aoo(fe) which can be used in evaluating Ad(b) using 
Bahcall & Soneira's (1980) procedure: 



A d (b) = A^b) 



exp 



— | d sin(6) | 
H 



(2) 



where b and d are the Galactic latitude and distance of the 
star, respectively. H is the scale height for the inters tellar 
dust which is adopted as 125 pc ((Marshall et al.ll2006h and 
Aoo(b) and Ad(b) are the total absorptions for the model and 
for the distance to the star, respectively. Then, the colour 



excess at the distance of the star, Ed{B 
ated using a specific form of Eq. 1: 

A d (b)=3AE d (B-V) 



V), can be evalu- 



(3) 



That value was used in Fiorucci & Munari's (2003) equations 
to obtain the total absorptions for the J, H and K s bands, 
i.e. Aj = 0.887 x E(B - V), A H = 0.565 x E(B - V) 
and Aks ~ 0.382 x E(B — V), which were used in Pogson's 
equation (mi — Ah = 51ogd — 5 + Ai\ i denotes a specific 
band) to evaluate distances. Contrary to the assumption 
above, the original (J — H) and (H — K s ) colour indices 
are not de-reddened. Hence, the application of the equations 
(1) to (3) is iterated until the distance d and colour index 
Ed(B — V) approach constant values. 

Our resulting distribution of E(B — V) reddening for 
our 21310 RAVE main sequence stars in the X-Y Galactic 
plane is given in Fig. 2. To analyze the reddening in the 
solar neighbourhood more accurately, we divided the stars 
into three subsamples according to their Galactic latitude: 
Fig. 2a shows the stars with |6| ^ 30°, in Fig. 2b the stars 
with 30° < |6| ^ 60° are shown, whereas Fig. 2c gives the 
stars with |6| > 60° . The first feature of the reddening distri- 
bution is the complex structure of the reddening in the first 
two panels. The local bubble, i.e. the region within 0.1 kpc 
distance from the Sun, is not affected by the reddening effect 
of the interstellar dust, whereas E(B — V) can be as high 
as 0.35 mag at larger distances. As expected high latitude 
stars (Fig. 2c) have smaller reddening values. The second 
feature is that one can not ignore the interstellar reddening 
even when using NIR bands. 

The distribution of distances (Fig. 3) shows that 80% of 
the sample stars have almost a normal distribution within 
the distance interval 0^d^0.4 kpc, whereas the overall dis- 
tribution which extends up to 1 kpc is skewed, with a median 
of 0.276 kpc. However, 97% of the sample stars are within 
d=0.6 kpc. 

The position of the sample stars in the rectangular co- 
ordinate system relative to the Sun is given in Fig. 4. The 
projected shapes both on the Galactic (X, Y) plane, and 
on the vertical (X, Z) plane of the sample show asymmetri- 
cal distributions. The median coordinates (X—60, F=-107, 
Z—-108 pc) of the sample stars confirm this appearance. The 
inhomogeneous structure is due to the incomplete observa- 
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Figure 3. Cumulative (a) and frequency (b) distributions of dis- 
tances of main sequence stars. 
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Figure 4. Space distributions of RAVE main sequence stars on 
two planes, (a) X-Y and (b) X-Z. 



tions of the RAVE project and that the programme stars 
were selected from the Southern Galactic Hemisphere. 



U, V and W are the components of a velocity vector of a 
star with respect to the Sun, where U is positive towards 
the Galactic centre (Z=0°, 6=0°), V is positive in the direc- 
tion of Galactic rotation (7=90°, 6=0°) and W is positive 
towards the North Galactic Pole (6=90°). 

Correction for differential Galactic rotation is necessary 
for accurate determination of U, V and W velocity compo- 
nents. The effect is proportional to the projection of the 
distance to the stars onto the Galactic plane, i.e. the W 
velocity component is not affected by Galactic differential 
rotatio n llMihalas fc Binnev 198ll ). We applied the proce- 



dure of lMihalas fe Binnev I l| l98ll ) to the distribution of the 



sample stars in the X-Y plane and estimated the first order 
Galactic differential rotation corrections for U and V veloc- 
ity components of the sample stars. The range of these cor- 
rections is -25.13<dt/<13.07 and -1.46«fV<2.28 km s" 1 for 
U and V , respectively. As expected, U is affected more than 
the V component. Also, the high values for the U compo- 
nent show that corrections for differential Galactic rotation 
can not be ignored. 

The uncertainty of the space velocity components 
U err , V err and W er r were computed by propagating the 
uncertainties of the proper motions, distances and ra- 
dial velocities, again using a (standard) algorithm by 
Ijohnson fc Soderblo m (1987). Then, the error for the total 
space motion of a star follows from the equation: 



(7err + Krr + Wen 



(4) 



The distributions of errors for both the total space mo- 
tion and that in each component are plotted in Fig. 5. The 
median and standard deviation for space velocity errors are 
Sen—8.01 km s _1 and s=16.99 km s _1 , respectively. We now 
remove the most discrepant data from the analysis, knowing 
that outliers in a survey such as this will preferentially in- 
clude stars which are systematically mis-analysed binaries, 
etc. Thus, we omit stars with errors that deviate by more 
than the sum of the standard error and the standard devi- 
ation, i.e. S err >25 km s" 1 . This removes 867 stars, 4% of 
the sample. Thus, our sample was reduced to 20453 stars, 
those with more robust space velocity components. After 
applying this constraint, the median values and the stan- 
dard deviations for the velocity components were reduced 
to (U err , V err , Werr)=(4.83±3.19, 4.40±2.82, 3.94±2.69) 
km s _1 . The two dimensional distribution of the velocity 
components for the reduced sample is given in Fig. 6. The 
most interesting feature for our present analysis in the (U, 
V) and (V, W) diagrams is the offsets of the zero points of 
the sample stars from the origins of the coordinate systems 
- these correspond to the solar velocity components. 



4 KINEMATICS 

We combined the distances estimated in Section 3 with 
RAVE kinematics and the available proper motions, apply- 
ing th e (standard) algorithms and t he transformation matri- 
ces o f lJohnson fc Soderbloml (|l987t ) to obtain their Galactic 
space velocity components (U, V, W). In the calculations, 
the epoch of J2000 was adopted as described in the Interna- 
tional Celestial Reference System (IC RS) of the Hipparcos 
and Tycho-2 Catalogues ilESAl 119971 ). The transformation 
matrices use the notation of a right handed system. Hence, 



5 POPULATION ANALYSIS 

We now wish to consider the population kinematics as 
a function of stellar population, using space motion as a 
statistical process to label stars as (probabilistic) "mem- 
be rs" of a stellar population. We u sed t he procedure 
of iBensbv. Feltzing fc Lundstroml (|2003l ) and iBensbv et al.l 
|2005h to allocate the main sequence sample (20453 stars) 
into populations and derived the solar space velocity compo- 
nents for the thin disc population to check the dependence o f 
LSR parameters on population. IBensbv et al.l (|2003l , I2005T ) 
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Figure 5. Error histograms for space velocity (a) and its compo- 
nents (b-d) for RAVE main sequence stars. The arrow in panel (a) 
indicates the upper limit of the total error adopted in this work. 
The shaded part of the histogram indicates the error for different 
velocity components of stars after removing the stars with large 
space velocity errors. 
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Figure 6. The distribution of velocity components of RAVE main 
sequence stars in two projections onto the Galactic Plane: (a) U- 
V and (b) W-V. 



assumed that the Galactic space velocities of stellar popu- 
lations with respect to the LSR have Gaussian distributions 
as follows: 



f(U, V, W) = k. exp 



UlSR {VlSR - Vasyr, 



Table 2. The space velocity component ranges for population 
types of RAVE main sequence stars. 



Parameters 


U 


V 


W 




(km s _1 ) 


(km s- 1 ) 


(km s _1 ) 


TD/D < 0.1 


(-85, 75) 


(-80, 40) 


(-55, 40) 


0.1 < TD/D s; 1 


(-105, 100) 


(-100, 50) 


(-70, 50) 


TD/D > 1 


(-165, 155) 


(-175, 60) 


(-115, 90) 



where 
k = 



(2^)3/2^ 



,<TV, 



><TWi 



(6) 



normalizes the expression. For consistency with other anal- 



yses, we adopt ay l 



and aw LS R as the character- 



istic velocity dispersions: 35, 20 and 16 km s for thin disc 
(£>); 67, 38 and 35 km s _1 for thick di sc (TP); 160, 90 an d 
90 km s _1 for halo (H), respectively jBensbv et alj|2003h . 
Vasym is the asymmetric drift: -15, -46 and -220 km s _1 for 
thin disc, thick disc and halo, respectively. Ulsr, Vlsr and 
Wlsr are LSR velocities. The space velocity components 
of the sample stars relative to the LSR were estimated by 
adding the values for the s pace velocity components evalu- 
ated bv lDehnen fc Binnevl |l99Sl) to the corresponding solar 
ones. 

The probability of a star of being "a member" of a given 
population is defined as the ratio of the f(U, V, W) distri- 
bution functions times the ratio of the local space densities 
for two populations. Thus, 

TD /P=^.?1£- T P/H=^.t^ (7) 

Xu Jd Xh Jh 

are the probabilities for a star of it being classified as 
a thick disc star relative to being a thin disc star, and 
relative to being a halo star, respectively. Xn, Xtd and 
Xh are the local space densities for thin disc, thick 
disc and halo, i.e. 0.94, 0.06, and 0.0015, respectively 



l|Robin et al"1ll996l: iBuser et all 19991 ). We followed the ar- 
gument of Bensbv et al.l ( 20051 ) and separated the sam- 



2al 

u LSR 



v LSR 



wj SR 

Wlsr 



pie stars into four categories: TD/D^O.l (high probability 
thin disc stars), 0.1<T_D/D^1 (low probability thin disc 
stars), 1<T_D/D^10 (low probability thick disc stars) and 
TP/P> 10 (high probability thick disc stars). It turned out 
that 18026 and 1552 stars of the sample were classified 
as high and low probability thin disc stars, respectively, 
whereas 412 and 463 stars are probabilistically thick disc 
and halo. The ranges and distributions of the galactic space 
velocity components for our different stellar population cat- 
egories are given in Table 2 and Fig. 7. 



6 RESULTS 

We analysed the U , V and W space velocity components 
for the main sequence sample (20453 stars) as well as for 
four subsamples, i.e. high probability thin disc stars of dif- 
ferent spectral types a) F, G and K (7V=17889 stars), b) F 
(A?=9654 stars), c) G (A r =5910 stars), d) K (iV=2325 stars) 
and estimated the modes of the Gaussian distributions for 
each category and each space velocity component. Spectral 
I J^S?} 38 were determined according to the 2MASS colours of 
the stars l|Bilir et alj|2008l . see Table 3) . The colour range is 
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Vfrms 1 ) V(kms') 

Figure 7. Distribution of main sequence stars of different popu- 
lation types on two Galactic space velocity planes: (a) U-V and 
(b) W-V. The contours indicate (outwards from the centre) high 
probability thin disc (TD/D ^ 0.1), low probability thin disc (0.1 
< TD/D ^ 1) and low and high probability thick disc and halo 
stars (TD/D > 1). 

0.08< J - #<0.30 for F spectral type and 0.30< J - #<0.42 
for G spectral type. The modes of these kinematic distribu- 
tion functions are the best estimators of the (negative) value 
of the solar motion relative to the LSR values for the stel- 
lar group in question. The histograms for the space velocity 
components for the main sequence star sample are given in 
Fig. 8. Gaussian fits, which are clearly an adequate descrip- 
tion of the data, reveal (7 =8.83±O.24, V = 14.19±O.34 and 
W =6.57±O.21 km s" 1 as the modes of the corresponding 
histograms. The modes of the histograms for the space ve- 
locity components for high probability thin disc stars of F, 
G and K spectral types, obtained by fitting them to the 
Gaussian distributions (Fig. 9), are a little different than 
the preceding set, i.e. (7 =8.5O±O.29, V = 13.38±O.43 and 
W Q =6.49±0.26 km s _1 . The largest difference in V can be 
confirmed by comparing Fig. 8b and Fig. 9b (see Section 7 
below for details). 

The distribution functions of Galactic space velocity 
components for high-probability thin disc stars, using our 
kinematic population assignment process, and presented 
separately for each of three spectral-type groups, are given 
in Table 3. We evaluated the modes for each distribution 
from the best-fit Gaussian. These results are also given in 
Table 3, together with those for the high-probability main 
sequence star sample and for the full sample. There are dif- 
ferences between the modes of the corresponding space ve- 
locity components for stars with F, G and K spectral types 
as well as between these modes and the corresponding ones 
for other sub-samples which will be discussed in Section 7. 



7 DISCUSSION 

By definition, the Local Standard of Rest (LSR) is the ref- 
erence frame of a hypothetical star at the galacto-centric 
distance of the Sun, and in the Galactic plane, that would 
be on a circular orbit in the gravitational potential of the 
Galaxy. We can attribute U , V and W velocity components 
of the Sun relative to this hypothetical star by observing the 
kinematic distribution function of stars at the solar neigh- 
bourhood, provided the stellar comparison sample we use is 
in dynamical equilibrium. We can similarly test the equilib- 



rium state of stellar populations by comparing their deduced 
solar peculiar motion, or equivalently their mean motion rel- 
ative to a LSR. One may wonder if the kinematically cold 
thin disc stars feel a large asymmetric gravitational poten- 
tial from, perhaps, an inner Galactic bar. Do spiral arms 
have systematic kinematic effects visible in young stars? Is 
the LSR defined by the high velocity halo stars, which feel 
largely extended dark matter potentials, the same as that 
of thick disc stars, and thin disc stars? These questions es- 
sentially repeat a query on the closeness of the higher-order 
Oort Constants to zero, for all Galactic components. 

The significance is such one wishes to determine the 
"solar peculiar velocity" for many samples, using several 
techniques. Doing this (see Table 1 above for a list) has 
led, in general to determinations of the peculiar motion ra- 
dially outwards from the Galactic centre Uq which are con- 
sistent and small, independent of the comparison sample 
chosen. Every stellar sub-population in the Galaxy seems 
remarkably radially stable and circularly symmetric, with 
no evidence for bar-like perturbations. Similarly, doing this 
for the net vertical motion, Wq provides the same conclu- 
sion. The sky fails to fall (see also Seabroke et al., 2008). 
The Galaxy, locally, seems not be significantly affected by 
an accreted stellar system, or asymmetric mass infall. The 
situation is less clear in the rotational direction, V . Here, 
there is less agreement between determinations, largely since 
model dependent assumptions need to be made to correct 
kinematics for energy balance, which transfers angular mo- 
mentum systematically into random motions with time. This 
may happen by slow diffusion, or by more dramatic radially- 
dependent effects. The relevant way forward is to determine 
Vq as a function of as many definable "populations" as is 
feasible, in many ways. We follow that approach here, where 
we find, in agreement with several recent studies, a value of 
Vq near 13 km s _1 . We determine this value with respect to 
three sub-groups of local main-sequence stars. These groups 
cover a very wide range of stellar ages. Yet they agree, sta- 
tistically, with each other, and with other determinations 
based on more local samples. The concept of a LSR again 
seems well-defined. 

We study stars which are thin disc dwarfs like the Sun. 
We expect a dispersion in the distribution of the velocity 
components due to their inherent distributions and differ- 
ences between the same velocity components of sample stars 
of different population types and ages. The way to reveal the 
corresponding velocity component is to fit the space motion 
data to a suitable distribution function. Encouragingly, a 
simple Gaussian provides an excellent approximation so we 
can identify the mean solar motion with its mode. This is 
the procedure we used in this work. 

The star sample was taken from RAVE, which extends 
much beyond the distance of the high-precision local Hip- 
parcos sample, i.e. 97% of the sample lies within the distance 
interval 0<d^0.6 kpc. We applied the following constraints 
to obtain a main sequence sample with relatively small er- 
rors: i) we selected stars with surface gravity 4<logg<5, ii) 
we omitted stars with (J — //)o<0.05 and (J — 7/)o>0.55 to 
avoid blue horizontal branch and red giant stars, iii) we ex- 
cluded stars with space velocity errors larger than 25 km s _1 , 
iv) we separated stars probabilistically according to their 
population types, and we used the high-probability thin disc 
sample to minimize contamination by thick disc/halo stars 



© RAS, MNRAS 000, [TV?? 



8 Co§kunoglu et al. 

Table 3. Space velocity components of the Sun with respect to the LSR for five population subsamples, as described in the text. 



Parameters 


Colour Range 


U 


V 


W 


N 








(km s — ) 


(km s — ) 


(km s — 




All sample 


0.05^(J - 


#) sS0.55 


8.83 ± 0.24 


14.19 ± 0.34 


6.57 ± 0.21 


20 453 


TD/D < 0.1 


0.05^(J- 


H) sS0.55 


8.50 ± 0.29 


13.38 ± 0.43 


6.49 ± 0.26 


18 026 


F Spectral Type 


0.08<(J - 


H) s:0.30 


8.35 ± 0.36 


13.14 ± 0.43 


6.24 ± 0.27 


9 654 


G Spectral Type 


0.30<(J - 


H) ^0A2 


9.25 ± 0.50 


14.42 ± 0.57 


6.67 ± 0.38 


5 910 


K Spectral Type 


0.42<(J - 


H) o ^0.55 


7.01 ± 0.67 


11.96 ± 0.66 


7.03 ± 0.38 


2 325 




-150 -100 -50 50 100 -100 -50 50 100 -100 -50 50 100 1 50 

U(kms- 1 ) V(kms- 1 ) W{km s" 1 ) 



Figure 8. The distribution functions of Galactic space velocities for the 20453 main sequence stars, with overlaid best-fit Gaussian 
distributions. 



as a preferred sample in our work. We estimated solar space 
velocity components for five samples: for all our RAVE main 
sequence stars; for the sub-sample of high-probability thin 
disc main sequence stars; and separately for the thin disc 
main sequence stars of F, G and K spectral types. The re- 
sults (see Table 3 above) are consistent, within sampling er- 
rors, across all sub-samples. They are also in agreement with 
recent determinations using other population subsamples, 
and over very different ranges. There is increasing agreement 
that the solar peculiar velocity in the rotation direction has 
been underestimated in older studies. The internal agree- 
ment between our different samples, and between this study 
and others, provides evidence that the Galactic potential 
near the Sun is symmetric, with no evident time dependent 
variation. 
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